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Abstract 



The calculation of the Kapitza boundary resistance between dissimilar har- 
"c^ ' monic solids has since long (Little [Can. J. Phys. 37, 334 (1959)]) suffered 

H . from a paradox: this resistance erroneously tends to a finite value in the limit 

— Jj . of identical solids. We resolve this paradox by calculating temperature dif- 

Ph ', ferences in the final heat-transporting state, rather than with respect to the 

j^ I initial state of local equilibrium. For a one-dimensional model we thus derive 

an exact, paradox-free formula for the boundary resistance. The analogy to 

ballistic electron transport is explained. 
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I. INTRODUCTION 

The occurrence of a thermal resistance at an interface between two materials, also 
known as Kapitzaa resistance, is a well-known phenomenon, in particular in low-temperature 
physics. In fact, the present work has been prompted in the course of an investigation of 
non-isothermal stochastic processes in Josephson junction devices.H The basic physical mech- 
anism is partial reflection of the phonons transporting the heat current. For the low tem- 
peratures (in comparison with the Debye temperature) at which the boundary contributes 
significantly to the total thermal resistance of a typical sample, phonon wavelengths are 
much larger than the lattice spacing and a continuum description of the materials should 
be adequate. In this case the reflection process becomes a standard problem in elasticity 
theory, and using semi-classical arguments to determine phonon densities and fluxes the 
Kapitza resistance was calculated along these lines in a classic paper by Little.H 

For an ideal interface between the two elastic media and an infinitesimal temperature 
difference between them so as to allow the use of linear response theory, a fully quantum 
statistical calculation of the transport process was given by Leung and Young.Q Besides 
confirming the intuitive expectation that for this idealized system the semi-classical result 
should be exact, the calculation is of interest in its own right as a non-perturbative evaluation 
of a Kubo formula. In addition, for a sample of finite size the dissipative behaviour implicit 
in the term 'resistance' was explicitly shown to hold only up to times at which Poincare 
recurrence sets in. 

In Refs. ^ and | the authors noted their result for the Kapitza resistance to be paradoxical 
in the limit of identical media: it tends to a finite value instead of vanishing. They attributed 
this artifact to the neglect of three-phonon scattering (which would take over as the dominant 
origin of resistance in this case), arguing that without such scattering the phonon mean free 
path is infinite so that it would become problematic to describe the system in terms of 
local temperatures. This explanation is unsatisfactory for the following reasons. First of 
all, three-phonon scattering is effective as a cause of heat resistance and thermalization 
only when Umklapp processes are possible.! Thus, one would have to endow the model 
with a lattice structure as well. However, as mentioned above, this should not be relevant. 
Secondly, the phonon mean free path is infinite for all values of the energy transmission 
coefficient T, and the authors of Refs. ^ and ^ rightly do not suspect their results in the case 
of small T. Thirdly, the vanishing of the boundary resistance for a homogeneous medium is 
a matter of principle, and should not depend on the relative magnitude of other resistance 
contributions. 

In this article we show that the aforementioned paradox is resolved in a simple manner 
by calculating temperature differences in the final, heat-transporting state of the system, 
and not with respect to the initial state of local equilibrium. For clarity of exposition we 
limit ourselves here to one space dimension. A report on the technically more involved 
three-dimensional case is currently in preparation.! In section |I| the corrected formula will 
be obtained by a semi-classical consideration in the spirit of Ref. |^. In Sec. |TT| this new 
result will be derived in the Kubo formalism. Apart from solving the Kapitza resistance 
paradox (for which Leung and Young's Kubo formula itself needs to be amended) the work 
of Ref. 1^ will be extended by giving the full space and time dependence of the heat current. 
In the one-dimensional case the relaxation to a stationary non-equilibrium state is thus 



shown to be degenerate, namely jump-wise. Also, one will see that the semi-classical result 
holds without the necessity of coarse-graining, partly because the problem has no length 
scale apart iroradhc (where c is the sound velocity!. The issue discussed in this paper has a 
striking analogyQa in ballistic electron transport .ciEl This analogy is the subject of Sec. |V. 
Some concluding remarks are made in Sec. |V]. 

II. SEMI-CLASSICAL ANALYSIS OF TEMPERATURE RELAXATION 



Consider two semi-infinite harmonic strings,aOO joined at x = 0, with sound velocities 
ci (cr) for X < {x > 0). For this system the boundary conductance K was defined in Ref. ^ 
as i^ = J'/ AT, where JT" is the energy current in the stationary state which is reached after 
the system has started at time t = in a state of local equilibrium with temperatures T^ 
for a; > and Ti = T, + AT for x < 0. Their result (Eq. (17)) is K = ^TT, with T 
given in ( |A2| ). The ensuing value for K in the case of identical strings (i.e., T = 1) is easily 
understood using a semi-classical picture (Ref. |^) in which energy- carrying fields propagate 
according to the wave equation. In view of the absence of an interface, excess energy at x < 
will then be radiated half to x = — oo, half to positive x. One thus arrives at K = ^cC, 
where C is the heat capacitance per unit length C = 7iT/3hc. This indeed gives the finite 
value K = ^T, in contrast with the expected result K = oo. The finite heat current is a 
consequence of the finite specific heat of the quantum string, combined with the finite speed 
of propagation (as opposed to diffusive transport). 

Now consider the general case of unequal strings, with initial thermal energy density e\ 
(et:) (= {'H\^r), where (H) is the expectation value of the Hamiltonian density; see Sec. |TTT|) 
for X < (x > 0) at time t = 0. At any t = ti > the situation is as in Fig. |l|. Properly 
accounting for left- and right-going radiation, the expression for the steady state value of 
the energy density e'^ {i = 1, r) is obtained by considering the history of an energy packet 
e[ 6xi and noticing that 5x^/5x\ = c^/c\. The result reads: 

, _ e\ + TZe\ + Te^Cr/c\ , _ Sy.^'Re^^ TeiC\/c^ 

with 7^ = 1 - r. Again putting Tx = T, + AT {AT < T^), e['^ = ttTP^ /Qhd —where we 
have used the specific heat as given above — , one gets 

T/ = Tr + (l-r/2)ATl 

t; = Tr + (T/2) AT j ^ AJ - /<, AJ . [Z) 

With the boundary conductance defined as the ratio of the stationary heat flux and the 
stationary temperature jump, one now finds 

K = ^ = ^^ = -T- (^) 

AT' AT AT' 6h n' ^ ^ 

which indeed diverges for T ^ 1. 



III. TEMPERATURE RELAXATION IN THE KUBO FORMALISM 



The result (|^) can also be calculated in the Kubo formalism. That calculation (Sec. [lllij| ), 
however, requires a nontrivial correction of the Kubo formula used in Ref. ^. In both 
Sees. Ill Al and IIII Bl we set h = 1. 



A. The Kubo formula 

At t = 0, to within first order in the temperature nonuniformity, the local equilibrium 
density matrix igi^ 



Plo 



Poll- p-^fdr ldx'6/3{^'){n{^', -ir) - {n{^ , -?r))o) I 



(4) 



where (3~^ is the equilibrium temperature, po = exp{—(3H) is the equilibrium density ma- 
trix {H is the full Hamiltonian) , 7i is again the Hamiltonian density and (■)q indicates an 
equilibrium average. Note that (H) is not restricted to the one- dimensional situation. When 
calculating the energy current (i7(x, t))^^^, the first and third term of pioc do not contribute, 
and one obtains 

{Ji^,t)),^^ = T-^j^ jdid5T{^){n{^, -iT)J{^,t)), (5) 

= (:r(x, 0))i„, + '- jd^ 5T(x') f^dt' Tr jpo f^dr e^-[7^(x'), i/]e-^V(x, t')\ . 

To begin with, (i7(x, 0)))^^^ vanishes because of time reversal invariance. In the second term, 
the trace can be taken as a single correlation function; one would then get (J7'(x, t))]^^ oc 
/q dt' {jQdTj'{x.', —iT)J'{x, t'))^ without an additional t' in the integrand. However, in order 
to relate {J^)^^^ to {JJ)^^^ we proceed differently and apply the Kubo identity: 



U(x,t))i„^ = ^|rfx'5r(x')/WTr{[po,H(x')]J(x,t')} 

/dx' 5T(x') Cdt' I ni^) + f rft"e"'^*"VJ^(x')e'''*", J^(x) 



I 
T 



(6) 



where the step taken on the second line of Eq. (j^) has been repeated. Integration by parts 
over x' and combining the t'- and ^''-integrals now yields 

(^(x, t))i„, = ^ /rfx' VT(x') J^dt' {t - t') ([J(x, t'), :r(x', 0)]), + 



it 
T 



rfx'5T(x')([7^(x'),J(x)])„ 



I /dx' VT(x') J^dt' {t - t') ( J(x, t'):r(x', 0))„,, + 



it 
f 



rfx'5T(xO([7^(x'),J(x)])„ 



(7) 



This Kubo formula will be crucial in Sec. [1IIB| . Apart from the boundary term, it differs 



from the one used m Ref. | m that the u.te'^d pjeseutly involves (' - * instead of ('.0 

To be sure, this modification does not affect the value of the stationary current. This 
follows from the explicit form for {JJ)odd calculated in the Appendix: 



{J{x, t)J{x\ 0))„,, = -{-^ r (t -a)--— 6'{t - a) + (a - -a) , (8) 




where x > 0, x' < 0, a = xjc^ — x' /c\, and where we returned to the one- dimensional case. 
Because the distributions 5"' and 5' are odd functions of their arguments, /q°° dt' {J{t')J)^^^ 
is seen to vanish. 

In the boundary term of the Kubo formula (|^), the commutator \H, J\ will be localized 
at x' = X. As will be shown explicitly in a simple case, this boundary term compensates 
a second unphysical term arising from the lower limit of the t'-integration in the first part 
of Eq. (^. In this first part, (j7'(x, t')j7'(x', 0))^^^ vanishes outside the 'sound cone'. Thus, 
for t' I only the point x' = x contributes in the x'-integration. In the calculation of 
the energy drop in Sec. |III B| , bT{x') vanishes identically near x' = x > and hence both 



unphysical contributions to the current are absent. To show their cancellation for arbitrary 
(5T(x') one needs the full space dependence of {JJ)^^^. Therefore, we will temporarily 
restrict ourselves to the case T = 1, where — by symmetry — formula (|^) does hold for any 
t and a. The equal-time commutator {[H^J]) can be calculated using Wick's theorem or, 
equivalently, by expanding the fields in creation and annihilation operators.li^ The answer 
is: 



-i\ 5"'{a) -Kd'ia) 
c \ 127r 3/?2 



mw,jm„--\"-^-^^\- (9) 



Inserting this in the second term of (0), and evaluating the first term of that formula to be 

5'{t — a) 71 ( ^, , ^, ,\ (5"{a) 7r(5(a)\l , , 

d{t-a)-6{-a)]-t\^^-—^\\ + {a^-a), 




247r 6/32 

(10) 

the unphysical terms which grow linearly in t are seen to cancel. 

B. The energy drop 

For the case of arbitrary T, the drop in energy density A(?-^)j^^ can be calculated from (|^ 
by using the detailed form of (H). However, only a few general properties of this correlation 
function will actually be needed (thereby suggesting the validity of the present analysis 
beyond the case of the ID elastic continuum). Namely, by inspection of the defining integrals 
it is easy to verify that for t, a > 0, {JJ)^^^ is an odd function /(s = t— a), which is localized 
at s = 0. As mentioned after Eq. (^, this immediately shows that the additional t-term 
in the Kubo relation (|^) gives a vanishing contribution to the current, except at t = a.li3 
On the other hand, omission of this term would lead to a A(7i)j^^ which is a factor two too 
large compared to the correct value obtained below in Eq. ([TT|). 



Let us now finally give a microscopic calculation of the change in energy density. Again 
choosing x > 0, one has 

J 

1 t°^ 

Cr Jo 

= lim ^^ rdt tdt' it - t') djit' - a) 

-r-oo Cr Jo Jo 

= lim ^— ^ / dt tf{-a) - / dt'f{t'-a) 

C-j. •^ (J L i' (J 

Const r , , . ,, , 

= Jim / dt{T- t)f{t - a) 

= iU(x,t = oo)),„,. (11) 

In the fourth and fifth line the first term vanishes because of the properties of / mentioned 
above. With {J)i^^ = |TAT, the value of A{TC)i^^ according to Eq. (|ll|) coincides with the 



one obtained in the semi-classical analysis of Sec. ||, leading to our result (^ 

IV. ANALOGY WITH BALLISTIC ELECTRON TRANSPORT 

The analysis of electric current in a one-dimensional conductor as a scattering prob- 
lem was pioneered by Landauer.ta Consider two well-separated reservoirs at electrochemical 
potentials fiy = ep + eV and /Xr = ep, respectively. Let the reservoirs be connected by a 
one-dimensional lead, interrupted by a barrier (at x = 0, say) with transmission coefficient 
T. Using the ID density of states l/hvp [vp being the Fermi velocity), at zero temperature 
the current per spin direction is obtained as / = T{evF){eV/hvF) = j^TV . This leads to a 

2 

conductance G = ^T in the case of a voltage measurement at the reservoirs. 

Near the barrier, however, the conductor is not well characterized by the reservoir poten- 
tials /ii^r- For energies between /xi and /Xr, at a position x < all left-going states are occupied, 
but only a fraction TZ of the right-going states is filled. This leads to an effective electro- 
chemical potential ii[ = Sy + eV{l + 7^)/2. Similarly, for x > one gets /ij, = ep + eVT /2. 
Calculating the conductivity using (/i( — /ir)/e = TZV as the relevant voltage difference for a 
measurement near the barrier, one finds the Landauer formula 

If one measures the voltage difference between the reservoirs, the finite upper limit of 
G = ^T as T I 1 makes perfect sense physically. The resulting G^^ = ^ has the significance 
of a purely geometric contact resistance between the reservoirs. On the other hand, G in 
Eq. (|12|) is the quantity which describes the intrinsic properties of the barrier. 

The analogy to the heat transport problem is now easily seen: the role of the electric 
current is played by the heat current, the electron densities from which the electrochemical 
potentials are determined are the analogue of the energy densities considered in Sees. || and 
m, and the effective electrochemical potentials translate into local temperatures. 

6 



The appropriate definition of ^[ is not trivial: one could argue that local chemical po- 
tentials are only defined up to the difference in potential pertaining to left- and right-going 
states separately. This leads to a conceptual uncertainty {AG)/G ~ 2T, and would there- 
fore allow a thermodynamic description of the current transport only if T -C 1. In that 
case the distinction between the two formulas for G obviously loses its meaning. The anal- 
ogous caveat holds for heat transport.EJ However, the potentials /ij^, leading to Eq. (|l^) are 
defined in such a way that, if the system would have been in local equilibrium at these 
potentials, the electron density would be the same as the actual density. This guarantees 
that the conductance ( ]T^ ) trivially satisfies the Einstein relation G = e^d^UpD, where the 
diffusivity D relates the particle current J^ to the particle density jump across the barrier 
Anp by J7p = -DAnp.B Moreover, this definition is the correct one in the following sense: 
if one couples two additional reservoirs infinitesimally weakly to the left and right of the 
barrier, the condition that no current will flow is that they must be at the potential /i( and 
/x^, respectively (see Ref. || for the precise statement). In our case one could — at least in 
principle — probe local temperatures by weakly coupling additional heat baths in a similar 
way. 

In the system described in Sees. || and |T|, the role of the reservoirs is played by the 



infinite leads themselves. For electron transport, this way of modeling the reservoirs was 
introduced earlier in Ref. |l^. The criterion for a good reservoir is that it should only re- 
emit electrons (phonons) after complete thermalization and loss of phase memory. The 
leads satisfy this condition by never re-emitting any particle. Violation of this condition, 
for example by taking finite leads, will lead to different results, as mentioned already in the 
Introduction. 

The absolute temperature, which enters as a prefactor in our formula for the boundary 
conductance, has no counterpart in the current transport problem. Its presence is a conse- 
quence of the Bose statistics obeyed by the phonons (see also the denominator in ( [Al|) ), as 
opposed to the Fermi statistics of electrons. 

V. CONCLUDING REMARKS 

The reasoning leading to the paradox-free formula (Eq. (^) for the Kapitza boundary 
conductance K calculated in this article is not invalidated by the presence of a thermalization 
mechanism like the weak three-phonon scattering contemplated in Refs. ^ and ^. On the 
contrary, such a mechanism rather brings the local energy distribution closer to a thermal 
one (under the condition that there should be a heat current J), and thus only strengthens 
the case for e[ as the energy from which stationary temperatures should be calculated. 

In conclusion: using the model of coupled harmonic strings, we have solved the Kapitza 
boundary resistance paradox (Little's formula) by means of an exact calculation of the 
properties of the interface in the final steady state, leading to the new formula Eq. (^. 
Besides, we have explained the analogy between local temperature measurements in the 
heat transport problem and voltage measurements in ballistic electron transport (Landauer's 
formula) . 
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APPENDIX: THE CURRENT-CURRENT CORRELATION FUNCTION 



We calculate 



With the Hamiltonian density Ti = ^tt^ + 



along the lines of Ref. 

^(V0)^ (where is the field operator, vr is its conjugate momentum, p is the mass density, 
and M is the modulus of elasticity), local energy conservation implies JT" = — |^{7r, V0} 
for the energy current^ Using this expression for JT, {JJ)^^^ can be expanded in terms 
of field-field correlation functions using Wick's theorem. By a KMS-relation, the field-field 
correlation functions can be obtained from the field-field commutators only, and these in 
turn follow from the classical wave equation. Carrying out this procedure leads to 



{J{x,t)J{x\Q)) 



odd 



T fdu rduj' ^'^ sin 



{uj + uj'){t-a) 



271 J 27r (l-e-^'^)(l-e-^'^') 
where T is defined in terms of the impedances Zi = piCi as 

AZiZ, 



+ [a ^ 



-a] 



(Al) 



r 



[z, + z,Y 



(A2) 



and X > 0, x' < as in the main text. The infrared divergence present in the field-field 
correlation functions because the energy of an oscillation vanishes in the long-wavelength 
limit, is seen to be absent in the current-current correlation function as such long-wavelength 
excitations hardly contribute to the energy transport. On the other hand, the ultraviolet 
divergence caused by taking the continuum limit for the string must be dealt with, e.g., 
by regularizing {J J) odd using a spatial smoothing. Taking a Lorentzian for that purpose 
leads to a cutoff e^^''^"'"'^ ' in the integral in ([Al]). However, since the integrand vanishes as 
r^ Q-P\^ I in ^j^e half plane lj + cj' < one may as well use the cutoff e~'^^^^^ •* . The integral 
then factorizes and the ensuing single integrals can be done exactly. With /3 = 1, this yields 



duj rdu' cu^'sin (cj + uj')t 



,-r]{uj+uj') 



2ti J 27r (l-e-'^)(l-e-'^') 



Im 



TT 



2sin\'n{r]-it)) 



{5*-2S5^}""FT» + '""[° (<' + "') 



17] 

^r(t)-^5(t) forr^iO. 



(A3) 



Scaling back to /5 7^ 1 and inserting this result into 
Eq. m. 



one obtains {JJ)r,dd ^^ given in 
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FIGURES 

FIG. 1. Two harmonic strings (with an interface at x = 0) at a time t = ti after the system 
started in a different local equilibrium on either side of a; = 0. 
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